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Abstract 

We suggest the diffuse approach to the relaxation processes within the kinetic theory for the 
Wigner distribution function. The diffusion and drift coefficients are evaluated taking into consid¬ 
eration the interparticle collisions on the distorted Fermi surface. Using the finite range interaction, 
we show that the momentum dependence of the diffuse coefficient D p (p ) has a maximum at Fermi 
momentum p = pf whereas the drift coefficient K p (p ) is negative and reaches a minimum at 
p ~ pf- For a cold Fermi system the diffusion coefficient takes the non-zero value which is caused 
by the relaxation on the distorted Fermi-surface at temperature T = 0. The numerical solution 
of the diffusion equation was performed for the particle-hole excitation in a nucleus with A = 16. 
The evaluated relaxation time r r ~ 8.3 • 10 -23 s is close to the corresponding result in a nuclear 
Fermi-liquid obtained within the kinetic theory. 

PACS numbers: 21.60.-n, 21.60.Ev, 24.30.Cz 
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I. INTRODUCTION 


The relaxation processes in many-body systems can be effectively studied within the 
kinetic theory, see Ref. [1] and references therein. The kinetic approaches operate with 
the kinetic equation, which is written for the distribution function in phase space. The 
advantage of the kinetic approaches is that the kinetic equation can be easily generalized 
to the case of finite temperatures. Certain difficulties arise when one tries to describe the 
relaxation and damping effects involving the collision integral in 9-dimension space 2-7]. 

To reduce the kinetic equation for the Wigner distribution function, we will follow the 
diffuse approach jsj] considering the relaxation on the distorted Fermi surface. The presence 
of the Fermi-surface distortion effects gives rise to some important consequences. Because of 
the Fermi-surface distortion, the scattering of particles leads to the relaxation and the damp¬ 
ing. The purpose of present paper is to study the diffusion and drift terms in a Fermi-system 
which can be applied to the description of the damping of collective and particle-hole exci¬ 
tations as well as the large amplitude dynamics. In Sec. 2 we consider the kinetic equation 
for the Wigner distribution function and reduce it applying the diffusion approximation. In 
Sec. 3 we establish the diffusion and drift coefficients taking into consideration the Fermi- 
surface distortion effects. The discussion of numerical results is presented in Sec. 4. Our 
conclusions are given in Sec. 5. 


II. COLLISION INTEGRAL WITHIN DIFFUSE APPROXIMATION 


We will restrict ourselves to the Born collision approximation in the kinetic equation for 
Wigner distribution function /(r,p;f) = / jg, 10]. Introducing the collision integral St{/}, 
we will write the kinetic equation as 


d/(r,pi;t) 

dt 


+ Lf(r,pi,t) = St {/}, 


( 1 ) 


where L is the driving operator. In the lowest order in h the driving operator L is given by 

Q 


L — — p • V r — CV r U) • V p . 

m 


( 2 ) 
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Here, the single-particle potential U includes, in general, the self-consistent and external 

fields. The collision integral St{/} in Eq. (0Q) can be written in the following form [ll] 

G . r f g 2 dp 2 dp 3 dp 4 (S u 
St{/} = ./ (2 «K)> W({Pi}) 

X /(Pl)/(P2)/(P3)/(P4) - /(Pl)/(P2)/(P3)/(P4) <5(Ae)<y(Ap) 

= j ffg, W^1(P1,P3)/(P1)/(P3) - ^ 3 (P1,P3)/(P1)/(P3) , (3) 

where g is the spin-isospin degeneracy factor, /(p) = /(r,p,t), /(p) = 1 — /(p), Ae = 
£i + e 2 — e 3 — £ 4 , Ap = pi + P 2 — P 3 — P 4 and W({pj}) is the probability of two-body 
collisions. The gain and loss terms W 3 tji(pi, P 3 ) in Eq. (J3]) are given by 

WW(Pi,Ps) = W (pi, P 3 ) = j W({pj})/(p 2 )/(p4)(5(Ae)(5(Ap), (4) 

Wb^PnPs) = ^(PnPs) = j w ({Pj})/(P2)/(P4)5(A£)(5(Ap). (5) 

The transition probability W p t^ q (p, q) in Eq. ((2D contains the square of the corresponding 
amplitude of scattering for the direct, 1 —» 3, and the reverse, 3 —)■ 1, transitions. The 
probability W p <=; q (p,q) includes also the distribution functions of the scattered particle in 
the initial and final states. We will assume that the main contribution to the scattering 
amplitude is given by the transitions which correspond to a small momentum transfer: 

|Pi — p 3 | -C Pfi where pp is the Fermi momentum, see also Ref. ll|. Introducing the new 
variables 

, 1 . . s 

s = P 3 -Pi and P = -( Pl + p 3 ) = Pl + -, 

we will apply the following expansions over small s : 

/(Ps) = /(Pi + s) ~ /(Pi) + s„V Pll „ /(Pi) + ^ s u s p W Pl ^W Pl:P /(pi), (6) 


W(p 1 ,p 3 ) = W(P,s) 


W (p 1; s) + i s u V pu „W (pi, s) + i s u s p V pitU V pup W (pi, s), 


and 


(7) 


W (pi, p 3 ) = IE(P. s) 


fR(Pi,s) + ^ s„X 7 pitU W (pi, s) + i s v s p S 7 Pl)V y pltP W (pi, s). (8) 
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Note that the spin-averaged probability of two-body collisions W({pj}) in Eq. (j3j) can be 
expressed in term of in-medium scattering cross section da/dD as 


W({Pi» 


2{2nh) 3 da 
m 2 dD 




(9) 


In the case of elastic collisions, the scattering cross section da/dD depends on the modulus 
square of the momentum transfer s only. 

Using the expansions Eqs. Q and (JHD ? one can reduce the kinetic equation (P) to the 
diffusion equation in the following form (see Eq. (1A9I) of Appendix A) 


dl 

dt 


Lf = -V, 


K p (p)/(p)/(p) +/ 2 (p)V p U) p (p) + V 2 [/(p)-Dp(p)], (10) 


where D p ( p) and K p ( p) represent the diffusion and drift terms, respectively. Both kinetic 
coefficients D p { p) and K p ( p) are derived by the following relations, see Appendix A, 

BvA p) = D p {p)8 v/t , D p ( p) = ^ J ^ s 2 W (p, s) (11) 

and 

K p {p)Vp,„e p = V PtP B vli ( p) - A v { p) = V PtU D p (p) - A v { p). ( 12 ) 


III. KINETIC COEFFICIENTS 


The obtained expressions (TTTh and (jT2T) for the kinetic coefficients imply a smallness of 
the momentum transfer s because of the expansions Eqs. ()6j) , (171) and (|8l) . To provide 
a small momentum transfer s in Eqs. (fill) and (1171) we will use the finite-radius inter¬ 
particle interaction with the following Gaussian form-factor v(r) = Vq exp(— r 2 /2r^) which is 
appropriate for calculations of the in-medium cross-section within the transport approaches 


13]. The differential cross section da/dfl in the first Born approximation is then given 
by (14| 

7rm 2 7> 0 2 


dD 


2ft 4 


■ exp (—4s 2 ro/ft 2 ) , 


(13) 


were ro and vq are the free parameters. 
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A. Diffusion term 


Using Eqs. (1A2j) and (TT3|) . we will rewrite the diffuse term D p { pi) of Eq. (TTT]) as 

2 6 2 r 

D p(P\) ~ J d P 2 dp 4 ds s 2 exp (-4s 2 r 2 //i 2 ) /(p 2 )/(p 4 ) 6 (p 2 - p 4 - s) 

x^£ 2 -£ 4 --). (14) 

V m / 

Integrating in Eq. ([Iff]) over s, we obtain (see Appendix B) 

„2 6 ? .2 C 00 A /' ^ 

D p (p ) « 247r2 0 / . 7 ° J dk k 5 J dx J dy exp [-8k 2 r%(l - xy)/h 2 ] 

x |(1 — xy) j ^8 k 2 rlVl — x 2 \/l — y 2 /ti 2 ^j — \Jl — x 2 \Jl — y 2 

x j (%k 2 rQ^/l — x 2 \Jl — y 2 /h 2 ^j | / ^ \Jk 2 + p 2 + 2 kpx^j 

X / ^\/fc 2 + p 2 + 2 kpy'j , (15) 


where x = cos 9 q and y = cos 6k- 

We will apply our consideration to the nuclear Fermi-liquid. For the numerical calcula¬ 
tions we will adopt the following parameters r 0 = 0.8 fm and Vq = —33 MeV, which provide 
a reasonable value for the in-medium nucleon-nucleon cross section cr to t — 20 mb. We will 
also use the Fermi distribution function 


f{p) = 


^1 + exp 


p 2 /2 m - A (T) 
T 


-1 


5 


(16) 


where T is the temperature, A (T) « e f [1 — (tt 2 /12){T/e F ) 2 ] is the chemical potential and 
£p = 37 MeV is the Fermi energy. 

The results of calculations of the diffusion coefficient D p (p) accordingly to Eq. (fT5]l with 
the Fermi distribution (1ff~6h are presented in Fig. 1. The calculations were performed for the 
different values of temperature T. As seen from Fig. 1, the momentum dependence of D p (p ) 
has the clearly observed maximum at the Fermi momentum p F for different temperatures. 
With an increase of temperature the diffusion coefficient increases as D p {jp = Pf,T) ~ T 2 . 
The temperature dependence of the diffusion coefficient D p (p = p Fl T) is shown in Fig. 2. 
As seen from Fig. 2, for a cold Fermi system the diffusion coefficient takes the non-zero value 
D(p F ,T = 0) = 2.1 • 1CU 22 MeV 2 -fm _2 -s which is caused by the relaxation on the distorted 
Fermi-surface in Eq. ([I]) which exists at T = 0 also. For high temperature regime where the 
Fermi statistic comes to Maxwell one, the diffusion coefficient D(p F ,T) behaves as a linear 
function of T in agreement with the Einstein’s fluctuation-dissipation theorem. 
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FIG. 1. Dependence of the diffusion coefficient D p (p ) on momentum p in units of pp, for temper¬ 
atures T = 0.1 MeV (solid line), 2 MeV (dashed line) and 4 MeV (dotted line). 



FIG. 2. Dependence of the diffusion, D p , and drift, K pi coefficients taken at p = pf on the 
temperature T. 

B. Drift term 

Using Eq. (jT2j) and e p = p 2 /2m, one can write the drift coefficient K p ( p) as 

A " p ( p ) m = ^p> vD p( p ) “ ^(P)’ ( 17 ) 
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FIG. 3. The same as in Fig. 1 but for the drift term I\ p (p). 


where 


A(p) = ~ [ 
P J 


gds 


PvSv W(p,s). 


(18) 


(2nh) 3 

For a spherically symmetric distribution f(p) the drift coefficient K p ( p) in Eq. (TT71) is 
reduced to the following form (see Appendix B) 

m / dD p (p ) 


k p(p) = 


A(p) 


(19) 


p \ dp 

where the first moment function A(p) is given by Eq. (IB7D of Appendix B. 

The presence of maximum of D p (p) at the Fermi momentum denotes that dD(p) / dp\ p=PF = 
0 and in agreement with Eq. (UUD the drift coefficient K p (p F ) is reduced as 


m 


K p (p F ) = - A(p F ). 

Pf 


( 20 ) 


In Fig. 3 we show the dependence of the drift coefficient K, p (p) on the momentum p in units 
of Fermi momentum p F . One can see that K p (p) has a minimum localized at p < p F which 
slightly depends on the temperature. Note that for p > p F the temperature dependence of 
K p (p ) is negligible. 

The numerical estimates of the diffuse coefficient D p (p) and the drift term K p (p) near the 
Fermi momentum give D p (p F ) ~ 10~ 22 MeV 2 -fm -2 • s and K p (p F ) ps —10” 23 MeV ■ fm -2 • s. 
Both obtained values of D p (p F ) and K p (p F ) agree with the phenomenological ones used 


earlier in Ref. 


15) 
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IV. NUMERICAL RESULTS 


Evaluating the transport coefficients D p (p ) and K p (p ) caused by the interparticle colli¬ 
sions on the distorted Fermi surface, we will restrict ourselves to a Fermi system which is 
homogeneous in r -space and omit the driving operator L in equation (ITQl) . Using the diffu¬ 
sion, D p (p), and drift, K p (p), coefficients from Eqs. (1T5T) and (fT9lh and solving the diffusion 
equation (11(1 with L = 0 one can evaluate the time evolution of the Wigner distribution 
function f(p,t). We will also assume a spherical Fermi surface of radius pp which is derived 
by the condition for the particle number A within a fixed volume V 


CPF 


47t gV 
(2n h) 3 


p 2 dp = A. 


The diffusion equation (fTOl) must be augmented by the initial condition for f(p,t)- We will 
consider the time evolution of the initial particle-hole (lplh) excitation which is derived at 


t — 0 as, see also Ref. 


id, 


finip) = [1 - 0(p - p[) + 9{p - p' 2 )] [1 - e(p - p F )] 

+ [1 -6{p-p 2 )\9{p-pi)9{p-p F ). (21) 


The distribution f\ n (p) of Eq. (j2TT) means the particle located at p\ < p < p 2 and the 
hole excitation at p\ < p < p' 2 for fixed pi > Pf and p 2 < pf , respectively. The intervals 
A p’ = p' 2 — p\ and A p = p 2 ~ Pi are derived from the conditions 



AngVdp 
(27 rh) 3 
AngVdp 
(277 h) 3 


P 2 fin(p,t = 0 ) 
P 2 fin(P,t = 0) 


A-l, 

1 . 


( 22 ) 


In Fig. 4 we have plotted the time evolution of the Wigner distribution function f(p,t ) 
for the initial particle-hole excitation in the Fermi-system with A = 16. We used the 
initial distribution / in (p, t = 0) given by Eq. (T2TT) and assumed the initial excitation energy 
E ex = 30 MeV. 

One can see from Fig. 4 that the momentum distribution f(p,t ) evolves to the Fermi 
equilibrium limit / eq (p) of Eq. (TT6l) with temperature T = \JE ex /a , where a = (tt 2 /6 )g(e F ) 
is the level density parameter 16] and g(ep) is the single-particle level density at the Fermi 
energy. The corresponding relaxation time r r for the diffusion process in Fig. 4 can be 









FIG. 4. Time evolution of initial distribution function (1211) in momentum space (in units of the 
Fermi momentum pf ) with p'i/pf — 0.71, p' 2 /pf — 0.75, Pi/pf — 1-12, P 2 /PF — 1-13, which 
corresponds to the initial excitation energy E ex = 30 MeV. The solid line is the equilibrium 
distribution / eq (p). 


obtained considering the time evolution of the deviation Sf(p,t) = f(p,t) — f eq (p ) of the 
distribution function f(p,t) from its equilibrium limit / eq (p). We introduce the mean square 
deviation 


A(f) = I dp [Sf(p,t)f 


dp df(p,t) 


= dp [Sf(p,t)f 


(23) 


The time dependence of A(t) for the distribution function f(p,t) from Fig. 4 is plotted in 
Fig. 5. The function A(f) can be fitted by the exponential dependence A (t) ~ exp(— t/r r ) 
which is shown in Fig. 5 as the dashed line. Both curves in Fig. 5 are normalized to the 
initial mean square deviation 

A 0 = J dp [5/ in (p)] 2 , 


where 5fi n (p ) = /; n (p) — fe q (p)- The corresponding relaxation time r r is estimated as T r 
8.3 • 10 -23 s. The obtained value of tv agrees with a typical collisional relaxation time t^cou 
in a nucleus r riCO n = 1CT 23 4- 10~ 22 s 


17 


19]. 
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FIG. 5. Dependence of the mean square deviation A (t) from Eq. (|23D on the time (solid line). The 
dashed line is the mean square fit to the exponential function A (t) ~ exp(—t/r r ). 


We note also that considering the time evolution of the distribution function to the 
equilibrium (solid line in Fig. 4), we reach the regime where the further change of distribution 
function is negligible. That means that the collision integral in equation (TTOl) should be zero 
in the case of equilibrium. To verify this fact we consider the time evolution of the collision 
integral (lA9j) . In Fig. 6 we show the ratio 


= I d P [St{/(t)}] 2 
f dp [St {f(t = 0)}] 2 ’ 


(24) 


which represents the time dependence of the mean square of collision integral St{/(f)} 
normalized to the initial value of St{/(t = 0)}. 

As seen from Fig. 6, the collision term St{/(£)} (1A9I) approaches to zero with time as it 
should be in the equilibrium limit 


V. CONCLUSIONS 

We have considered the kinetic approach with the collision integral for the Wigner dis¬ 
tribution function using the diffusion approximation. We have demonstrated that the colli- 
sional Landau-Vlasov kinetic equation can be reduced to the form of the diffusion equation. 
We have established the expressions for both the diffusion coefficient D p (p) and the drift 
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FIG. 6. Time dependence of the ratio Rst(t). 

term K p (p) for the Fermi-systems where the relaxation processes occur on the distorted 
Fermi surface. 

We have found out that the simplest isotropic assumption for the collision probability in 
the collision integral is insufficient for the correct description of the kinetic coefficients. The 
forward scattering should be ensured for the inter-particle interaction. For the description 
of the inter-particle interaction we have used the finite-range potential. Such kind of inter¬ 
particle interaction provides the necessary smallness of the momentum transfer and thereby 
the forward scattering of particles. 

For the nuclear Fermi-liquid we have calculated dependencies of the diffuse coefficient 
D p (p) and the drift term K p (p) on the momentum p and the temperature T. For the diffuse 
coefficient D p {p) the momentum dependence possesses the clearly observed maximum of 
order of 10~ 22 MeV 2 -fm _2 -s near the Fermi momentum for different temperatures. With an 
increase of temperature the value of the diffusion coefficient increases also. Excluding the 
small momentum values, the drift term K p (p) has negative sign and the minimum of order 
of 1(T 23 MeV • fm _2 -s. Inside the Fermi sphere for p < p F the drift term K p (p) shows the 
strong dependence on the temperature. For p > p F the temperature dependence of K p (p ) is 
practically negligible. 

We have established the finite value of the diffusion coefficient in a cold Fermi system with 
D(pp, T = 0) = 2.1 • 10~ 22 MeV 2 • frn -2 • s which is caused by the relaxation on the distorted 
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Fermi-surface at T = 0. The diffusion process was investigated numerically assuming the 
initial non-equilibrium Ipl/i-excitation for a finite Fermi-system with number of particles 
A = 16. It was shown that the lpl/i-excitation relaxes to the equilibrium Fermi distribution. 
The numerical estimate for the relaxation time r r gives r r ~ 8.3 • 10 _23 s which is close to 
the corresponding estimates in a nuclear Fermi-liquid obtained within the kinetic theory. 


Appendix A: Moment representation for the collision integral. 


Applying the expansions of IF(pi,Ps) and !F(pi,P 3 ) over small s = p 3 — pi (see Eqs. 
o and (ED), the collision integral d3J) can be written as 

stum ! gds 


(2nh) 3 ' W ( Pl> S ) “ W ( pi > s ) ) /(Pi)/(Pi) 

Su { (iT (Pi, s)/(pi) + IT (pi, s)/(pi)) V Pl> „/(pi) 
^/(Pi)/(Pi)V p1i , (^IT(pi, s) - IT(pi,s) 

|lT(pi, s)/(p 1 )V pi ^V pi)A 1 /(pi) 

/( Pl )Vp li i/IT(p 1 , s) + f( Pl )V pi ,„W( Pl , s)) V pl)/ J( PlJ 
/(Pi) (w ( P i, s)V Pli j/V PliAt / ( Pl ) 


i /(Pi)V pi ,,V pi)/i (IT( Pl , s) - IT( Pl , s) 


Using then Eqs. 


and ([HD, we obtain 

d P2 d P J{ P2 )f( P 4)S (p 2 - P4 - s) 

P-s' 


IT(P, s ) ^ ^^(s 2 ) 

m z d\ l 


and 


x 5 ( £ 2 — £4 


2 g da / _ 2 


m 


X 5 ( £ 2 - £4 


d P2 d P4 /(p 2 )/(p4)5 (P 2 - P4 - s) 

P s 


m 


The collision integral of Eq. (lAljl is transformed to the following form 
gds 


St{/} = 


(27 Th)' c 


x V Pl)i ,IT(pi, s)} + 


{lT(pi, s) (/(pi) - /(pi)) V Pl ,,/(pi) + /(pi)/(pi) 
SvS p |lT(pi, s)/(pi) V Pl)i ,V Pl , M /(pi 


gds 
(27r/yd 


P 1,V 


(IT( Pl , s)) V pliM /( Pl ) + IT( Pl , s)/( Pl )V pliI/ V pliAt /( Pl )} 


(Al) 


(A2) 


(A3) 


(A4) 
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Introducing the two first moments A u (pi) and -B w/ 1 (pi) of scattering probability W( pi, s) 


as 

4,(pi) = j s » W(Pi>s)> ^m(Pi) = \f { l % 3 s » s n w (pi,s), 

the collision integral (1A4|) is reduced as 


St{/} = V P1)V 


^(pi)/(pi)/(pi) + ^( Pl )V piiA J(pi) 


Using the relations 


(AS) 


(A 6 ) 


^(pi)V Pl , M /(pi) 

= V pii/1 (^p(pi)/(pi)/(pi)) - /(pi)/(pi)V pl># 1 S I/p (p 1 ) + ^(p 1 )V pi , M / 2 ( Pl ) 

= V pi ,p (B I/ p(pi)/(pi)) - / 2 (pi)V Pl)At i?j //i (pi) - /(pi)/(pi)Vp ljM -B^(pi), (A7) 


as 


we rewrite the collision integral 

St {/} = - V pi)i , (Vp^iUP!) “ ^(Pi)} /(Pi)Z(Pi) + / 2 (Pi)V pliMJ B i// ,(p 1 ) 
A V PliV V PliAt (/(pi)-B^(pi))■ 


(A 8 ) 


Finally, we obtain 

St{/} = 


PI,*' 


Ap(pi)(Vp 1 ,^ pi )/(pi)/(pi) + / 2 (pi) V pii j / -Dp(p 1 ) 


V 2 


Pl,V 


[/(Pi)-Dp(Pi)] 


(A9) 


Appendix B: Diffuse and drift terms 

Integrating in Eq. (fT4l) over s, one obtains 

A>(Pt) ~ 9 ^Sh 7 J dp ' 2cip4 ' P2 _ Pl ) 2exp ( 4 (P 2 - P4) 2 Vo/ft 2 ) /(p 2 )/(p 4 ) 

x 6 fe 2 - £4 - — (p 2 - p 4 )) • (Bl) 

The argument of the delta function in Eq. (IBID can be written as 

£ 2 - £4 - — (P2 - P4) = ((P2 - Pi) 2 - (P4 - Pi) 2 ) , 

m 2m 

where £j = p 2 /2 m. We will introduce the new variables P 2 — Pi = q and p 4 — pi = k which 
allows one to rewrite Eq. (IBID in the following form 

D P {Pi) « 9 <£^2 fit J dq dk ( q _ k ) 2 exp ( _4 ( q “ k ) 2 r o/h 2 ) /(q + Pi)/(k + Pi) 

x S (q 2 — k 2 ) . (B 2 ) 
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Using the spherically symmetric distribution function /(p) and the relation 

5(x — a) + 5{x + a) 


5 ( x 2 — a 2 ) = 


2\a\ 


we 


will rewrite Eq. (1B2D 


as 


kdk / (q 2 + k 2 — 2 qk cos qk 


x exp (—4(g 2 + A; 2 — 2gfc cos qk)rg/h 2 ) {h(g — Ac) + S(q + fc)} 


x / \/V + Pi + 2qpi cos qpi 


)/(^ 


k 2 + pf + 2kpi cos kpi I . 


(B3) 


Integrating in Eq. 

D p {pi) ' 


over q, we obtain 


g 2 mrQVQ 


247 T 2 h 7 


dk k 5 / dUfc / dfL (1 — cos qk 


'o 

exp (^—8k 2 rl(l — cosqk)//i 2 j / 
f k 2 + p\ + 2kpi cosk-pi^ . 


x exp ( —8A; rg(l — cos qk )/h ) / k 2 + pf + 2kpi cos q pi 


x 


Using the addition theorem for spherical 
coordinate system 


(B4) 

harmonics, we can write in an arbitrary spherical 


Taking 
and using 


cos qk = cos 6 q cos 6k + sin 6 q sin 6k cos (cp q + c fa ), 
cos qfpi = cos 6 q cos 6 P1 + sin 6 q sin 6 pi cos((p q + (j ) pi ), 
cos k-pi = cos 6k cos 6 Pl + sin 6k sin 6 Pl cos (cftk + (ft Pl )• 

into account the spherically symmetry of the distribution functions /(pj) = f(Pj) 

cosk-px = cos 6k, 


coscfpi = cos 6 q , 


we obtain 


D p {px) « 9 2 ™2ff ° J q dk k5 f d ^ k f d ~^ ( X “ cos q k ) 

x exp 8fc 2 ro(l — cosqk)/h 2 j / k 2 + p\ + 2kpi cos 6 q 
x / [^Jk 2 + p\ + 2kp 1 cos 6k^j . 


(B5) 


We will also introduce the following notations 


p2n p2tt 

j(a) = d(j) q d4>k exp (a(4> q + 4 > k )), 
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p2n /*2n 

j{a) = d(p q d(j) k cos( 0 g + 0 fc ) exp [a(<j> q + (p k )), 

Jo Jo 

which are attributed as functions of the argument a. Using j(a) and j(a), one can write 
the angular integrals in Eq. ( 1 B 5 D as 


f 2 tt 


r2ir 


d:0q 


d(p k (1 — cos qk ] exp ( — 8 Arrg(l — cos qk)/ h 


= exp (— 8k 1 2 rl(l — cos 8 q cos 9 k )/K 2 ) 

x | (1 — cos 0 q cos 0 k )j(8k 2 r-g sin 9 q sin 9 k / h 2 ) — sin 9 q sin 9 k j(8k 2 ry sin 9 q sin 9 k /H 2 ) 

Finally, we obtain the diffuse coefficient D p (p ) in the form given by Eq. (1131) . 

The drift coefficient K p ( p) of Eq. (1T2|) can be reduced similarly to above described 
procedure. Using derivation of Eq. (fT8]h we will reduce A(pi) to the following form 

7t(pi) ~ J dq dk Pl (q - k) exp (-4(q - k) 2 r 2 /fr 2 ) /( q + Pi)/(k + p x ) 

X 5 (q 2 - k 2 ) , 


or 


A(p i 


g 2 mrlvl f°° 

87 T 2 H 7 Jn 


dk k 4 / dVt k / dkl q (cos9 q — cos9 k ) 


x exp ( — 8k 2 rl(l — cos qk)//i 2 ) / ( \ k 2 + p\ + 2 kpi cos 9, 


x / (^Jk 2 + p 2 + 2kpi cos . 


Finally, we obtain 
A(p) 


g 2 mrQVQ 


8ir 2 n 7 j dk kA f dx J d V ( x “ y) 


x exp (— 8 fc 2 ro(l — xy)/h 2 ) j (8k 2 r^Vl — x 2 \Jl — y 2 /hr 
x / (^\/k 2 + p 2 + 2 kpi'X^ f (y\/k 2 + p 2 + 2kpy S j . 


(B 6 ) 


(B 7 ) 
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